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1.0  INTRODUCTION 

The  calculation  of  two-dimensional,  turbulent,  separated  flow  using  the  boundary-layer 
equations  has  achieved  some  degree  of  success,  as  reported  in  Refs.  1  through  5.  Provided 
solutions  of  the  boundary-layer  equations  yield  acceptable  engineering  accuracy,  their  use  to 
compute  flows  with  separation  is  attractive  because  of  simplicity  and  cost  reduction 
compared  to  solution  of  the  full  Navier-Stokes  equations.  Boundary-layer  computation 
methods  can  be  classed  as  either  integral  or  differential  techniques.  In  general,  integral 
methods  are  simpler  and  require  fewer  computational  resources  than  do  differential 
methods  because  of  the  built-in  empiricism  such  as  a  velocity  profile  representation  and 
auxiliary  relations  for  shape  factors,  skin  friction,  etc.  This  report  presents  a  means  of 
analytically  describing  separated  turbulent  boundary-layer  velocity  profiles.  The  velocity 
profile  representation  is  used  to  develop  a  skin  friction  and  shape  factor  correlation  for 
turbulent  separated  flow. 

2.0  ANALYTICAL  DESCRIPTION  OF  SEPARATED 
BOUNDARY-LAYER  VELOCITY  PROFILES 

The  equation  proposed  to  describe  reversed  flow  velocity  profiles  on  smooth, 
impermeable,  adiabatic  walls  is  an  extension  of  the  expression  presented  by  Whitfield  (Ref, 
6),  which  is  a  composite  function  of  the  form 


which  consists  of  an  inner  expression  originally  presented  in  Ref.  7  and  an  outer  expression 
derived  in  Ref.  6.  The  inner  solution  as  presented  in  Ref.  6  is  given  by 

u+  =  — ! —  tan'l  (0.09  y!+'l 

'  009  (2) 


where 


(3) 


_+ 

V 


ury 

v 


(4) 
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ur  =  (Cf/2)/j  %  (5) 

and  the  subscript  “i”  refers  to  the  inner  region.  Assumptions  made  in  deriving  Eq.  (2)  for 
attached  flow  were  as  follows;  (l)  the  total  shear  stress  (molecular  plus  turbulent)  is  constant 
and  equal  to  the  wall  value,  ?w;  (2)  the  density  is  constant  in  the  wall  region;  (3)  the  turbulent 
kinetic  energy  is  proportional  to  the  Reynolds  stress,  -  <u 'v'  >;  and  (4)  the  turbulent  kinetic 
energy  (or  -  <  u 1  v '  > )  is  an  explicit  function  of  u  + .  These  assumptions  are  also  made  herein 
to  retain  the  form  of  Eq.  (2). 

The  slope  of  the  velocity  profile  at  the  wall,  computed  from  Eq.  (2),  is  positive  for  Cr  > 
0.  A  negative  slope  for  separated  flow  can  be  obtained  by  taking  the  velocity  distribution  as 

u  +  =  — —  tan'1  (0.09  y+)  (6) 

1  0.09  V  ' 

and  defining  ur  as 

»r  =  (|Cf|/2)‘4ire  (7) 


The  outer  expression  is  derived  such  that  the  complete  solution,  Eq.  (1),  approaches  the 
correct  asymptotes  at  the  wall  and  infinity.  From  Eq.  (6),  as  y  —  oo,  u,+  —  -7r/0. 18. 
Therefore,  for  u+  to  have  the  correct  limiting  value  of  u*as  y  —  oo,  u+  must  behave  as 


(8) 


as  y  —  oo.  In  addition,  because  ufbehaves  correctly  for  small  values  of  y,  u+  must  behave 
as  u+  -  0  as  y  -  0.  Therefore,  in  terms  of  the  outer  variable,  y/0,  the  form  considered  in 
Ref.  6, 


- 

_ + 

/  n  \ 

fy  \ 

uo  = 

ue  - 

\  0.18/ 

gi 

i 

is  retained,  where  g(y/0)  behaves  as  g(0)  -  0  and  g(oo)  =  l. 


(9) 


In  Ref.  6,  the  form  of  g(y/0)  was  determined  from  the  analytical  relation  for  attached 
flows,  as  follows: 
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(10) 


This  relationship  was  plotted  versus  y/0  for  numerous  experimentally  measured  velocity 
profiles.  The  resulting  curves  were  analytically  fit  by  the  relation 


(ID 


where  a  and  b  are  parameters  that  are  functions  of  Cf,  H,  and  Reg.  This  same  functional 
form  (Eq.  11)  is  retained  in  the  present  work.  Therefore,  the  analytical  representation  of 
two-dimensional,  separated,  turbulent  boundary-layer  velocity  profiles  becomes 


U 


-----  tan'l  (0.09  v+)  - 
0.09 


-+ 

u  - 

(-—) 

,  >/ 
tanh  '2 

o 

V  o  18/_ 

IAi 

(12) 


where 


-+  _  Refl  .  I. 


-  _+ 

U  U 


The  only  difference  between  the  expression  derived  in  Ref.  6  (for  attached  flow)  and  Eq.  (12) 
is  the  appearance  of  a  negative  sign  in  the  coefficients  of  the  trigonometric  functions  which 
is  a  consequence  of  modifying  Eq.  (2)  so  that  a  negative  velocity  slope  occurs  at  the  wall  for 
separated  flow.  Thus,  a  convenient  method  to  extend  the  result  obtained  in  Ref.  6  for 
attached  boundary-layer  velocity  profiles  to  include  profiles  with  reversed  flow  is  to  define 


S  = 


c 


f 


(13) 


and  rewrite  Eq.  (12)  as 
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_+ 

u 


— tan'1  (0.09  y+)  +  / u  + - —  \  tatih^ 

0.09  \  c  0.18  / 


(14) 


where  the  parameters  a  and  b  are  functions  of  Cf,  H,  and  Reg. 

3.0  PROFILE  CORRELATIONS 
3.1  FITTING  OF  VELOCITY  PROFILES 

Experimentally  measured  velocity  profiles  were  used  to  determine  the  parameters  a,  b, 
and  u+  which  appear  in  Eq.  (14).  The  determination  of  u+ provides  a  means  of  estimating 
skin  friction  for  separated  flows  in  a  manner  similar  to  that  used  by  Clauser  (Ref.  8)  for 
attached,  incompressible,  turbulent  boundary  layers.  Unfortunately,  only  two  experiments 
involving  detailed  flow  measurements  in  a  separated  turbulent  boundary  layer  could  be 
found  in  the  literature:  the  experiment  of  Simpson  et  al.  (Ref.  9),  which  involved  an 
incompressible  flow  with  an  airfoil-type  pressure  distribution,  and  the  experimental 
investigation  of  Alber  et  al.  (Ref.  10),  who  probed  the  transonic  flow  over  a  bump  on  a  wind 
tunnel  wall. 


Rearranging  the  argument  of  the  inverse  tangent  term  of  Eq.  (14), 


_+ 

y 


d 


(15) 


yields 


where 


(16) 


(17) 


Also,  the  argument  of  the  hyperbolic  tangent  function  of  Eq.  (14)  can  be  written  as 


(18) 
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where 


a 


(19) 


Therefore,  Eq.  (14)  becomes 


tan 


0.09 


1  T  /  Re\__ 

0  09  (it  )y 

L  \"e/  - 


s*  \ 
ms ) 


lanh 


'■  [V] 


(20) 


The  important  consequence  of  the  latter  manipulations  is  that  the  velocity  ratio  u/ue  (u/ue 
=  u+ /u|)  is  now  a  function  of  the  physical  distance,  y,  local  unit  Reynolds  number,  and 
local  skin  friction  (incompressible).  Because  the  available  experimental  data  do  not  have 
significant  compressibility  effects,  the  approximation  Re  •  y  *  Re  •  y  was  used.  Therefore, 
given  u/ue  at  three  points  in  the  boundary  layer,  a  system  of  three  equations  [Eq.  (20)  for 
each  point]  with  three  unknowns  [skin  friction  (in  terms  of  u~),  a,  and  b]  is  obtained.  With 
u+  a,  and  b  known,  one  can  obtain  d  using  Eq.  (20)  and  determine  a  from  Eq.  (19). 

The  experimentally  measured  velocity  profiles  of  Simpson  (Ref.  9)  and  Alber  (Ref.  10) 
were  fit  using  the  above  procedure;  that  is,  for  each  measured  profile,  Eq.  (20)  was  solved 
iteratively  (by  Newton’s  method)  for  u+  a,  and  b.  These  parameters  are  different  for 
each  velocity  profile.  Achieving  the  best  fit  of  a  given  velocity  profile  was  a  trial  and  error 
process  which  consisted  of  changing  the  three  points  needed  to  fit  the  profile.  Note  that  u+ 
is,  by  definition,  always  positive.  For  the  separated  profiles,  a  negative  skin  friction  was 
inferred  by  knowing,  a  priori,  that  the  measured  velocity  near  the  wall  was  negative. 
Therefore,  for  fitting  the  velocity  profile,  Eq.  (13)  was  modified  to  be 


Comparisons  of  fitted  and  measured  velocity  profiles  (separated  and  attached)  for  both 
experimental  investigations  are  presented  in  Fig.  1.  The  agreement  is  reasonably  good 
throughout  the  entire  range  of  y.  This  was  the  case  for  nearly  all  the  experimentally 
measured  profiles  considered  for  evaluation.  For  each  separated  profile,  the  best  fit  resulted 
from  using  the  maximum  negative  velocity  ratio  as  the  first  point  from  the  wall. 
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Upstream  of  separation 
(Sta.  -9.0) 


b.  Within  separation  (Sta.  4.5) 


c.  Downstream  of  separation  (Sta.  7.0) 


Figure  1.  Comparison  of  measured  attached  and  separated 
boundary-layer  velocity  profiles  with  Eq.  (14). 


d.  Upstream  of  separation  e.  Within  separation  (x-139.1  in.) 
(x=60  in.) 


f.  Within  separation  (x=157.1) 


Figure  1.  Concluded. 
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It  should  be  noted  that  the  measurements  of  Alber  (Ref.  10)  were  in  the  transonic  regime 
whereas  Eq.  (20)  was  developed  for  incompressible  flows.  However,  it  was  reported  by 
Lewis,  Kubota,  and  Webb  (Ref.  11)  that  for  compressible  flows  (up  to  Me  =  8.18),  the 
relation 


(22) 


is  suitable  for  relating  compressible  variables  to  their  incompressible  counterparts. 
Furthermore,  Whitfield  (Ref.  6)  shows  that  to  a  reasonable  approximation 


where  0k  refers  to  the  kinematic  boundary-layer  momentum  thickness  defined  by 


(24) 


Therefore,  Eq.  (20)  is  approximate  for  compressible  flows  if  rewritten  as 


where 


and 


-  + 


> 


(26) 


(27) 


and  Re9  which  appears  in  Eq.  (26)  is  obtained  by  applying  Coles’  “Law  of  Corresponding 
Stations”  (Ref.  12), 

cf  Hefl  =  c:(  Refl  (28) 
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in  conjunction  with  the  correlation  offered  by  Winter  and  Gaudet  (Ref.  13),  which  relates  Cf 
to  Cf  by  the  relation  Cf  =  Cf  [1  +  (M2/5)]l/2. 


3.2  SKIN  FRICTION 


It  is  of  interest  to  compare  skin  frictions  as  determined  from  Eq.  (20)  to  those  reported  in 
Refs.  9  and  10.  This  comparison  is  presented  in  Fig.  2,  which  shows  skin  friction  versus 
incompressible  shape  factor,  H.  Simpson  (Ref.  9)  measured  Q  directly  by  flush-mounted 
hot  film  devices,  whereas  the  skin  friction  reported  by  Alber  (Ref.  10)  was  inferred  from  a 
least-squares  fit  of  the  measured  velocity  profile  data  to  Coles’  law  of  the  wall  and  law  of  the 
wake  (see  Ref.  14).  The  skin  friction  data  reported  in  Ref.  10  have  been  transformed  to 


2.0r- 


1.5 


1.0 

'’S 

%  0.5 
io 


-0.5 


Alber  etal.  (Ref.  10)  Data 


l 


I  nferred  from  Fit  of  Eq.  (20)  to 
Experimental  Data 
I  nferred  from  Fit  of  Coles'  Laws  of 
Wall  and  Wake  to  Experimental  Data 


Simpson  (Ref.  9)  Data  j 

Solid  Line  Represents  Eq.  30)  for  Ree  =  21,500. 


I  nferred  from  Fit  of  Eq.  (20)  to 
Experimental  Data 

Measured  (Flush-Mounted  Hot  Film) 


>  Zw'°o  «bB 

o  * - CO - - 


-1. 


1 

2 


1 

5 


6 


7 


Figure  2.  Correlation  of  incompressible  skin  friction, 
Cf,  for  separated,  incompressible  flow. 


equivalent  incompressible  values  as  suggested  by  Winter  and  Gaudet  (Ref.  13)  according  to 
the  relation  Cf  =  Cf  [1  +  <M2/5)]1/2.  Simpson  (Ref.  9)  reported  that  precise  interpretation 
of  his  measured  results  was  difficult  downstream  of  intermittent  separation  because  the 
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flush-mounted  hot  film  probes  used  were  directionally  insensitive.  The  differences  between 
skin  frictions  computed  by  the  present  method  and  those  reported  in  Ref.  10  are  attributed 
to  the  fact  that  the  basic  velocity  profile  representations  are  different;  i.e.,  Eq.  (20)  or  Coles’ 
composite  form  of  the  law  of  the  wall  and  law  of  the  wake.  The  shape  factor,  H,  used 
in  Fig.  2  was  obtained  by  numerically  integrating  (using  Simpson  integration)  the 
profile  fit  used  to  obtain  skin  friction  fromjhe  simultaneous  solution  of  Eq.  (20)  at  three 
points  across  the  layer.  The  shape  factor,  i-T,  is  defined  by 

— - - -  (29) 


For  the  data  obtained  in  compressible  flow  (Ref.  10),  recourse  was  made  to  the 
approximation  u/ue  =  u/ue,  y/6  —  y/d k,  and  thus  =  H. 

The  solid  line  in  Fig.  2  represents  the  analytic  expression  derived  to  approximate  the 
numerical  and  experimental  results  for  Re,  =  21,500  [the  average  momentum  thickness 
Reynolds  number  of  the  Alber  experiment  (Ref.  10)].  The  expression  selected  was 


C 


f  = 


-  1 


(30) 


The  form  of  the  second  term  of  Eq.  (30)  was  chosen  such  that  for  shape  factors  less  than 
about  2.3,  the  expression  degenerates  into  the  relationship  given  by  White  (Ref.  15).  Any 
dependence  of  the  incompressible  skin  friction  on  Reynolds  number  has  been  assumed  to  be 
accounted  for  by  the  first  term  of  Eq.  (30)  [for  H  >  3.5,  Re,  has  little  effect  on  Cf, 
according  to  Eq.  (30)]. 


3.3  SHAPE  FACTOR  CORRELATION 

The  parameters  a  and  b  which  appear  in  Eq.  (14)  are  determined  by  matching  the  velocity 
distribution  at  two  points,  y/$  =  2  and  5  (Ref.  6  discusses  why  y/6  =  2  and  5  were  chosen 
for  the  match  points).  These  parameters  were  correlated  with  H  (Ref.  6)  by  plotting  u/ue  at 
y/6  =  2  and  5  versus  H  for  numerous  measured  velocity  profiles  (Ref.  14).  However,  the 
correlations  established  in  Ref.  6  did  not  include  separated  flow.  The  behavior  of  u/u£  with 
H  at  y/6  -  2  and  5  is  presented  in  Fig.  3  for  both  attached  and  separated  flow.  The  data 
were  compiled  in  Ref.  6  from  the  Stanford  conference  (Ref.  14)  and  inferred  from  fitting  the 
velocity  profile,  Eq.  (14),  to  the  measured  velocity  data  of  Simpson  (Ref.  9)  and  Alber  (Ref. 
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10).  The  dashed  lines  represent  the  correlations  established  in  Ref.  6,  and  the  solid  lines  are 
the  present  results  derived  by  including  the  data  inferred  from  the  separated  profiles.  The 
correlation  of  u/ue  at  y/0  =  2, 


U 

e 


1 

1 .95 


tanh"1 

(85  H\  0.364 

V  7.5  / 

(31) 


represents  the  large  shape  factor  data  reasonably  well  (Fig.  3a)  whereas  the  proposed 
variation  of  u/ue  with  H  at  y/0  =  5, 

^  (5)  =  0.155  +  0.795  sech  ^  0.51  (R  -  1.95) J  ^ 

falls  below  the  data  inferred  from  the  Alber  experiment  (Fig.  3b)  (Ref.  10)  for  H  >  3  but 
represents  the  data  inferred  from  the  Simpson  experiment  (Ref.  9)  reasonably  well.  Both 
relationships  approximate  the  Whitfield  correlations  (Ref.  6)  for  H  <  3.  It  should  be  noted 


a.  y/0  =  2 

Figure  3.  Correlation  of  experimental  velocity 
profile  data. 
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b.  y  10  =  5 
Figure  3.  Concluded. 


that  the  B  used  in  determining  the  velocity  ratio  of  y /&  =  2  and  5  was  obtained  by 
numerically  integrating  the  fitted  “raw”  profile  (dimensional  y  versus  u/ue)  using  the 
relation 


6  = 


(33) 


as  opposed  to  using  the  experimentally  determined  momentum  thickness.  This  was  done  so 
that  the  resulting  correlations  would  be  consistent  with  the  velocity  profile,  Eq.  (14).  The 
relationships  between  compressible  and  incompressible  variables,  Eqs.  (22)  and  (23),  were 
used  when  the  compressible  data  of  Alber  (Ref.  10)  were  correlated. 


Equations  (31)  and  (32)  can  be  used  to  generate  separated,  incompressible,  turbulent 
boundary-layer  velocity  profiles  from  H,  0+  and  Ree.  The  procedure  is  outlined  in  Table  1. 

As  derived  in  Ref.  16,  there  are  four  integral  length  scales  which  appear  in  the  boundary- 
layer  momentum  and  mean  flow  kinetic  energy  integral  equations: 


(34) 
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and 


In  addition,  the  following  shape  factors  were  defined  in  Ref.  16: 


(35) 

(36) 


(37) 


(38) 


H 


e* 


§1 

6 


(39) 


S 

e 


(40) 


these  shape  factors  were  also  correlated  as  a  function  of  H  and  boundary-layer  edge  Mach 
number,  Me.  For  the  present  work,  attention  will  be  restricted  to  establishing  a  relationship 
between  H3.  and  H  for  incompressible  flow,  Me  *  0.  As  a  consequence,  the  correlations 
derived  in  Ref.  16  for  H6.  =  H*.  (H,  Me)  and  H^.  =  Hs«  (H,  Me)  remain  unchanged 
because  for  Me  ~  0,  H$+  is  approximately  equal  to  H  and  H&*.  is  approximately  equal  to  0. 


The  velocity  profile,  Eq.  (14),  was  used  to  carry  out  the  numerical  integrations  (using 
Simpson’s  integration)  necessary  to  establish  the  correlation  between  H#*  and  H.  Because  Cr 
=  Cf  (H,  Re#),  Eq.  (30),  the  velocity  profile  is  only  a  function  of  H  and  Re@.  As  pointed  out 
by  Whitfield  (Ref.  16),  the  influence  of  Re*  on  the  velocity  profile  is  small  compared  to  that 
of  H .  The  results  presented  are  for  a  representative  turbulent  momentum  thickness  Reynolds 
number,  Re9  =  50,000. 

Figure  4  presents  the  numerical  results  for  Hp*  =  H@*(H)  as  open  symbols.  The  analytic 
expression 

=  1.48061  +  3.83781  e""^  r  0.33 - —  Lan'1 

0  8  5184 


iu"-M  -  t 

1  23 


-^0.33  -  — lanM  { 1.2874  x  10"6)  (lO7'11  ) 


1  4  57  6  I 


(41) 


was  chosen  to  represent  the  numerical  results. 
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Table  1.  Summary  of  Procedure  for  Computation  of  Separated 
Turbulent  Boundary-Layer  Velocity  Distributions 


lt-ci{LiH<Miienl 


I'q,  IS 
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Open  Symbols  Represent  the  Numerical  Results  tor 
Eq.  1141  Using  the  Velocity  Profile  Given  in  Table  1 


Figure  4.  Correlation  of  Hg*  for  adiabatic, 
incompressible  flow. 

4.0  CONCLUDING  REMARKS 

An  analytical  expression  that  reasonably  describes  the  velocity  distribution  in  a 
separated,  turbulent  boundary  layer  has  been  developed.  The  expression  is  an  extension  of 
the  form  derived  by  Whitfield  (Ref.  6)  for  attached  flows  on  smooth,  impermeable, 
adiabatic  walls.  The  analytical  expression  was  fit  to  available  experimentally  measured, 
separated  velocity  profiles  to  establish  velocity  correlations  for  u/ue  at  y/v  =  2  and  5  and  to 
establish  a  local  skin  friction  correlation  in  the  form  Cf  =  Cf  (H,  Rev),  which  allows  Cr  to 
become  negative.  These  correlations  and  the  analytical  velocity  profile  expression  were  used 
to  obtain  a  reversed  flow  boundary-layer  shape  factor  correlation.  The  analytical  velocity 
profile  expression  depends  on  local  values  of  skin  friction,  shape  factor,  and  Reynolds 
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number  based  on  momentum  thickness,  and  describes  attached  or  separated  flow  over  the 
entire  y  domain  (0  <  y  <  oo). 

Because  of  the  relatively  small  amount  of  experimental  data  which  were  available,  the 
relations  developed  cannot  yet  be  considered  universally  applicable  to  all  flow  situations. 
The  correlations  should  be  improved  as  more  experimental  data  become  available. 
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NOMENCLATURE 

Parameter  in  Eqs.  (11),  (12),  (14),  (18),  and  (19) 

Defined  by  Eq.  (19) 

Parameter  in  Eqs.  (11),  (12),  (14),  (18),  and  (19) 

Local  skin  friction  coefficient,  2rw/eeu^ 

Function  defined  by  Eq.  (11) 

Incompressible  shape  factor,  defined  by  Eq.  (29) 

Shape  factor  based  on  h*,  defined  by  Eq.  (38) 


a 

a 

b 

C, 

g 

H 

H6. 
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Shape  factor  based  on  5**,  defined  by  Eq.  (40) 

He*  Shape  factor  based  on  8 *,  defined  by  Eq.  (39) 

M  Mach  number 

Re  Local  unit  Reynolds  number,  ue/ 

Ree  Local  momentum  thickness  Reynolds  number,  u <&/vt 
S  Parameter  defined  by  Eq.  (13) 

u  Mean  velocity  in  the  axial  direction 

u+  Boundary-layer  velocity  coordinate,  u/uT 

uj4-  Inner  solution  for  u  + 

u+  Outer  solution  for  u  + 

uT  Friction  velocity,  (|Cf|/2)l/2  ue 

x  Coordinate  along  body  surface 

y  Coordinate  normal  to  body  surface 

y+  Boundary-layer  y  coordinate,  ury/v 

6*  Boundary-layer  displacement  thickness,  defined  by  Eq.  (34) 

6**  Boundary-layer  density  thickness,  defined  by  Eq.  (37) 

8  Boundary-layer  momentum  thickness,  defined  by  Eq.  (35) 

8  Equivalent  incompressible  momentum  thickness,  defined  by  Eq.  (33) 

8*  Boundary-layer  energy  thickness,  defined  by  Eq.  (36) 

0k  Boundary-layer  kinematic  momentum  thickness,  defined  by  Eq.  (24) 

v  Kinematic  viscosity 

e  Density 

7  Total  shear  stress 
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SUBSCRIPTS 

e  Boundary-layer  edge  value 

i  Inner  region  of  a  boundary  layer 

o  Outer  region  of  a  boundary  layer 

w  Wall  value 

oo  Infinity  or  free-stream  value 

SUPERCRIPTS 

—  Denotes  low-speed  or  incompressible  value 

'  Denotes  fluctuating  quantity 

SPECIAL  NOTATION 

<  >  Indicates  time-averaged  quantity,  e.g.,  <□%,>  =  —  f'  ( nv)  dL ' 

L-»ou  2t 
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